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The Maxwell equations for the spherical components of the electromagnetic fields outside sources 
do not separate into equations for each component alone. We show, however, that general solutions 
can be obtained by separation of variables in the case of azimuthal symmetry. Boundary condi- 
tions are easier to apply to these solutions, and their forms highlight the similarities and differences 
between the electric and magnetic cases in both time-independent and time-dependent situations. 
Instructive examples of direct calculation of electric and magnetic fields from localized charge and 
current distributions are presented. 
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Las ecuaciones de Maxwell para las componentes esfericas de campos electromagneticos en regiones 
libres de fuentes no son separables en ecuaciones para cada una de sus componentes. Se muestra, 
sin embargo, que soluciones generales pueden ser obtenidas por separation de variables en el caso 
de simetrfa azimutal. Las condiciones de borde son faciles de aplicar para estas soluciones, y sus 
formas destacan las similitudes y diferencias entre los casos electrico y magnetico, tanto para las 
situaciones independientes del tiempo como para las de tiempo dependientes. Se presentan ejemplos 
instructivos de calculos directos de campos electricos y magneticos producidos por distribuciones 
localizadas de cargas y corrientes. 
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1. INTRODUCTION 

The Maxwell equations for the electromagnetic field 
vectors, expressed in the International System of Units 
(SI), are Q 

V-D = ( o, VxE (>L> 



dt 



V ■ B = 0, 



VxH = J - 



<9D 
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where the source terms p and J describe the densities 
of electric charge and current, respectively. For a linear, 
isotropic medium D and H are connected with the basic 
fields E and B by the constitutive relations 



D = eE, H = B//x, 



(2) 



where e and p, are the permittivity and permeability of 
the medium, respectively. 

The boundary conditions for fields at a boundary sur- 
face between two different media are [2| 



n ■ (Di 



Do 



Ps, 



nx(Ei -E 2 ) = 0, 



drawn from the second into the first region. The interior 
and exterior fields satisfy the homogeneous vector wave 
equations 



V 2 E 



V 2 B 



efj, 
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which are obtained from Eqs. (JIJ and J5J) for regions 
free of charge and current by combining the two curl 
equations and making use of the divergence equations 
together with the vector identity 



V 2 ( ) = V(V- ) - Vx(Vx 



(5) 



Changes in the electromagnetic fields propagate with 
speed v = \jJep, . 

Without any loss of generality, we may consider only 
harmonic time dependence for sources and fields: 



p{r,t) = p(r)e~ 



E(r,i) = E(r)e- 



J(r,t) =J(r)e' 



B(r,i) = B(r)e" 



(6) 
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where ps and Jg denote the surface charge and current 
densities, respectively, and the normal unit vector n is 



where the real part of each expression is implied. Equa- 
tion (0J then becomes time- independent: 



V 2 E + fc 2 E = 0, V 2 B + fc 2 B = 0, 



(7) 
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where k 2 = efiui 2 . These are vector Helmholtz equations 
for transverse fields having zero divergence. Their solu- 
tions subject to arbitrary boundary conditions are con- 
sidered more complicated than those of the correspond- 
ing scalar equations, since only in Cartesian coordinates 
the Laplacian of a vector field is the vector sum of the 
Laplacian of its separated components. For spherical co- 
ordinates, as for any other curvilinear coordinate system, 
we are faced with a highly complicated set of three si- 
multaneous equations, each equation involving all three 
components of the vector field. This complication is well 
known and general techniques for solving these equa- 
tions have been developed, based on a dyadic Green's 
function which transforms the boundary conditions and 
source densities into the vector solution Q. We shall 
show, however, that in the case of spherical boundary 
conditions with azimuthal symmetry, the solution can be 
obtained more conveniently by means of separation of 
variables. Several applications of physical interest can 
then be treated in this simplified way, so avoiding the 
dyadic method 0. 

Actually, the usual technique for solving boundary- 
value problems introduces the electromagnetic potentials 
as intermediary field quantities. These are defined by Q 



B = VxA. 



E 



-V6- 



dA 



(8) 



The purpose of this paper is to get general solutions 
of the electromagnetic vector equations in spherical co- 
ordinates with azimuthal symmetry using separation of 
variables in spite of having equations that mix field com- 
ponents. Boundary conditions are easier to apply to 
these solutions, and their forms highlight the similari- 
ties and differences between the electric and magnetic 
cases in both time-independent and time-dependent sit- 
uations. The approach shows that boundary- value prob- 
lems can be solved for the electric and magnetic vector 
fields directly, and that the process involves the same 
kind of mathematics as the usual approach of solving for 
potentials. The material in this work may be used in a 
beginning graduate course in classical electromagnetism 
or mathematical methods for physicists. It is organized 
as follows. In Sec. 2 we describe the method for the 
static case showing how the mathematical complications 
of solving the vector field equations are easily overcome 
by means of separation of variables. In Sec. 3 we ex- 
tend the method to discuss the case of time- varying fields. 
Concluding remarks are given in Sec. 4. 



2. STATIC FIELDS 

For steady-state electric and magnetic phenomena, the 
fields outside sources satisfy the vector Laplace equations 



with the subsidiary Lorentz condition 



V 2 E = 0, V 2 B = 0, 



(13) 



V • A + -£ = 0. 

at 



(9) 



It is then found that these potentials satisfy the inhomo- 
geneous wave equations 

d 2 4> p 



(10) 



which together with the Lorentz condition form a set 
of equations equivalent to the Maxwell equations. The 
boundary conditions for the potentials may be deduced 
from Eq. JSJ. 

For fields that vary with an angular frequency to, i.e. 



(r,i) = 6{v)e 



A(r,i) = A(r)e" 



fill 



we get equations that do not depend on time in regions 
free of charge and current: 



k 2 6 = 0, 



W 2 A + k 2 A = 0, 



(12) 



where only transverse components with zero divergence 
are involved. Supposing all the charge and current are on 
the bounding surfaces, solutions in different regions can 
be connected through the boundary conditions indicated 
in Eq. ©. To demonstrate the features of the treat- 
ment, we first consider boundary-value problems with 
azimuthal symmetry in electrostatics. The solution of 
stationary current problems in magnetostatics is mathe- 
matically identical. 

Combining the expressions for Vx(VxE) = and 
V • E = in spherical coordinates and assuming no in- 
dependence, we find using Eq. that the components 
of the electric field E r and Eg satisfy the equations 



(V 2 E) r = 



(V 2 E) 



1 d 2 
~ 2 

i d dE r 



r 2 sin 8 89 
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= 0, (14) 
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Equation i|14[l is for E r alone, whereas Eq. (|15fl involves 
both components. There is also a separated equation for 

E,„: 



which are like those in Eq. J7J for the electric and mag- 
netic induction fields, so that in general we also confront, 
for the vector potential, the mathematical complexities 
mentioned above for the electromagnetic fields. 



(V 2 E), = \ ^(rE,) 



1 



d 



r 2 sin 86 
1 



r 2 sin 2 ( 



E v = 0. (16) 



3 



In this paper, however, we will not be concerned about 
those cylindrical symmetry cases where only the tp- 
component of the vector field is nonzero because a scalar 
technique of separation of variables is already known to 
obtain the solution 6]. 

Using the transverse condition 

V • E = 1 J-(r 2 P r ) + iUin 9 Eg) = 0, (17) 

r z or r sin oO 

where azimuthal symmetry is assumed, Eq. ()14|l implies 

(18) 



d ( w \ dEr n 



which is consistent with Eq. (|15fl . Thus, to obtain Eg 
from E r we can consider either Eq. Ijl7|l or Eq. (|18J) . 
These equations correspond to choosing a gauge when 
this method is applied to the vector potential. 

Now, in order to solve Eq. Ijl4|l for E r , we refer to the 
method of separation of variables and write the product 
form 



£;(r,9) = ^P(fl), 



(19) 



which leads to the following separated differential equa- 
tions: 



d 2 u n(n + 1) 
dr 2 r 2 



u = 0, 



(20) 



-^^ sin ^ +n ( n+1 ) P = ' ( 21 > 
sin dO at) 

where n(n + 1) is the separation constant. The solution 
of Eq. flSE) is 



u{r)=ar n+1 + — , 



(22) 



where a and b are arbitrary constants. Equation 121|) 
is the Legendre equation of order n and the only solu- 
tion which is single valued, finite and continuous over the 
whole interval corresponds to the Legendre polynomial 
P n (cos 6), n being restricted to positive integer values. 
Thus the general solution for E r is 

E r (r, 0) = f^ (dnr*- 1 + ) P«(cos 9). (23) 

The simplest way of solving Eq. H17|) for Eg is to use the 
series expansion 

Eg(r, 6) = J2 Mr) 3^«(cos 6), (24) 



n=0 



where v n (r) are functions to be determined. By replacing 
Eqs. (gSJ and JUJl into Eq. |T7|), it is found that 



v n {r) = — r 
n 



+ 1 r™+ 2 



(25) 



for n > 1 with a a — 0; this null factor in Eq. (|23|l means 
the absence of static field terms of the 1/r type, which 
are in reality typical of radiative fields as shown below. 
Clearly, the solutions given in Eqs. (J23J), (|23} and 
satisfy Eq. I|18|) . The coefficients a n and b n are to be de- 
termined from the boundary conditions. For complete- 
ness, we include here the well-behaved general solution 
of Eq. ijlUj) : 



EJr,( 



X 

E 

71=0 



P„(cos 9). (26) 



Thus, Eqs. 123|) - (|26|l formally give all three components 
of the electric field. The same type of equations applies 
in magnetostatics. However, the boundary conditions of 
Eq. J2J will make the difference, implying in particular 
that b Q = in the series expansion of Eq. I|23|) in magne- 
tostatics; this being primarily related to the absence of 
magnetic monopoles. 

To illustrate the use of the above formulas, we consider 
the simple example of the electric field due to a ring of 
radius a with total charge Q uniformly distributed and 
lying in the x-y plane. It is usually solved through the 
scalar potential method by using the result of the po- 
tential along the z-axis obtained from Coulomb's law [3 • 
The surface charge density on r = a, localized at 9 = ir/2, 
is written as 



Ps{a,9) 



Q 

2-ixa 2 



<5(cos 9), 



(27) 



which may be expanded using the well-known Legendre 
series 



<5(cos 9) = 



2n+ 1 



n=0 



P„(0) P„(cos 9), (28) 



with P n (0) given by 



P 2n+1 (0) = 0, P 2 „(0) 



(-l)"(2n+l)! 



2 2n (n 
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(29) 



Taking into account the cylindrical symmetry of the sys- 
tem, and the requirement that the series solutions in 
Eqs. I|23() - (|25(l have to be finite at the origin, vanish at 
infinity and satisfy the boundary conditions of Eq. © 
at r = a for all values of the angle 6, namely, Eg con- 
tinuous at r = a and E r discontinuous at r — a, it is 
straightforwardly found that the spherical components 
of the electric field are 



E r ir,9) = 



Q 



J2 P n(0) P«(C0S 



n=0 



(. + !)(!)' 



r > a 



"HJ ' r<a 



(30) 
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E g (r,6) = - 
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4ire r 2 



£>„(0)pi(cos 9) 
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£p^(0)P^(cos 
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r > a 



n+l 



(31) 



r < a 



and J5y, = 0, where P^(cos 9) = [d/dff) P„(cos 0) is an 
associated Legendre function. Note in particular that the 
coefficient b in Eq. becomes Q/Att€ for r > a, as 
expected. Also, the discontinuity of the nth component 
of E r in Eq. H3()|l at r = a is connected according to 
Eq. (J3J with the corresponding component of the surface 
charge density ps obtained from Eqs. I|27|) and (|28p. ex- 
hibiting the unity of the multipole expansions of fields 
and sources (see Ref. ||). 

To clarify the application of the formulas in the case of 
magnetostatics and also compare with electrostatics, we 
consider next the magnetic analog of the above example, 
that is, the magnetic induction field from a circular cur- 
rent loop of radius a lying in the x-y plane and carrying a 
constant current /. The surface current density on r = a 
can be written as 



Js(a, 6, ip) = — S(cos 9) <p, 
a 



(32) 



where for the delta function is now convenient to use the 
expansion 



n + 1 \r 

1 / r xn+2 



r > a 



r < a 



(36) 



Note that, as anticipated for magnetostatic problems, the 
coefficient b in Eq. (|23|l is equal to zero. Also, as ex- 
pected, the discontinuity of the nth component of Bg in 
Eq. i|36|) at r = a is connected according to Eq. JJJ with 
the corresponding component of the surface current den- 
sity Js v obtained from Eqs. (|32[) and H33[). Another char- 
acteristic difference with the electrostatic analog is that 
the coefficient P„(0) appears instead of P n (0). This can 
be traced to the fact that the inhomogeneous boundary 
condition, as given by Eq. J3J, is applied to the angular 
component of the magnetic induction field in Eqs. 1241) - 
(|25|l . as opposed to the corresponding inhomogeneous 
boundary condition acting on the radial component of 
the electric field in Eq. {23$ . The fields in Eqs. 
(|36|l are usually obtained through the vector potential 
method by using the expression of the magnetic induc- 
tion field along the z-axis calculated from the Biot and 
Savart law 0. An alternative technique is mere integra- 
tion of the vector potential 0- Our treatment has the 
advantage of introducing a considerable simplification on 
the procedure of applying the boundary conditions on 
the magnetic induction field directly. 



<S(cos 6) = ^T 



2n + 1 



-< 2n(n+ 1) 

n— v 7 



P^(0) Pi (cos 9), (33) 



which follows from the completeness relation for the 
spherical harmonics after multiplication by e~ t(p and in- 
tegration over if. The values for P^(0) are 

PL(0) = 0, P 2 V 1 (0)= ( ~ 1) 2 t n 1 ( y ( 34 ) 

Because of the cylindrical symmetry of the system, B v — 
0. By requiring that the field be finite at the origin, 
vanish at infinity and satisfy the boundary conditions of 
Eq. at r = a, the series solutions in Eqs. I|23|l - H25|) 
for the magnetic case lead to the following radial and 
angular components of the magnetic induction field: 



B r (r,9) = - 



2r 3 



2 00 



5>»P n (cos 



x < 



n=0 



n+2 



r > a 



r < a 



(35) 



TIME- VARYING FIELDS 



By using Eqs. (TJ, and © it is seen that outside 
sources the fields are related by 



(37) 



so that we only need to solve Eq. JJJ for B. Alternatively, 
we can solve for E, and obtain B through the expression 



B = VxE. 



(38) 



In the following, we choose to deal with the Hclmholtz 
equation for the magnetic induction field. The reason is 
to exhibit similarities and differences with the static case 
treated in Sec. 2. 

In the case of spherical boundary surfaces with az- 
imuthal symmetry, the B r and Bg components of the 
magnetic induction satisfy the following equations: 



(V 2 B) r + k 2 B r 
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d 

d9 {sm 



(r 2 B r ) 
„ dB r 



1 



09 



sin 
- k 2 B r 



0, (39) 
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, 2 N 2 1 a 2 , x 1 d 2 B r 

(^B)g + k 2 Bg = - 

+ k 2 Ba = 0. 



(40) 



Similarly, for the B v component we would have the equa- 
tion 



(V 2 B) v + fc 2 B v = \ -^{rB v ) 



1 



d 



r 2 sin 9 89 



x (sm 
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2- B v + k 2 B v = 0. (41) 



These are analogous to Eqs. <|14l) . I|15|l and 116() in con- 
nection with the vector Laplace equation. In order to 
solve Eq. we let 



B r (r,9) = 3 -^-P(6), 
r 



whereupon separation yields 



d 2 3 , 2 dj 



dr 2 



dr 



i{n + 1) 



3 = 0, 



(42) 



(43) 



and Eq. I|21|l . where the constant n(n + 1) is the sepa- 
ration parameter. Equation i)43|l is the spherical Bessel 
equation of order n with variable kr. Therefore, the gen- 
eral solution for B r is 



B r (r,t 



E 



jn{kr) 



n n (kr) 



P n (cos9). (44) 



Depending on boundary conditions, the spherical Hankcl 
functions hn ' instead of the spherical Bessel functions 
j n , n„ may be used. For Bg we again write 



B g (r,0) = ^Tw n (r) — P„(cos 



71=0 



(45) 



and use V • B = to obtain now 



— [r j n (kr)] 



n(n + l)r dr 

— [rn n (kf)], 



n(n + l)r dr 



(46) 



for n > 1 with a = 6 = 0. The other coefficients a n 
and 6„ are determined so that the boundary conditions 
for the vector field are exactly satisfied. In the case of 
the B v component, the general solution is 

oc ^ 

B v (r,9) = J2l 



n=0 



The same type of equations applies for the electric field. 

As an example, we shall consider the problem of the 
magnetic induction field from a current I = I e~ lut in 
a circular loop of radius a lying in the x-y plane. It is 



the time- varying version of the case solved in Sec. 2. The 
surface density current on r = a is then 

J s (o, 9, <f,t) = — S(cos 9) e~ wt 0, (48) 
a 

which can be expanded using Eq. l|33|l . The complete se- 
ries solution of the Helmholtz equation for the magnetic 
induction field, which is finite at the origin, represents 
outgoing waves at infinity and satisfies the boundary con- 
ditions of Eq. @ at r = a, becomes 



B r (r,9,t) = -i 



Holoka 



(2ti + 1)P„ 1 (0) 



j n (ka) h[l\kr) 
xP„(cos 9) { (49) 
Jn (kr) hiP(ka) 



B g {r,9,t) = 



xP^(cos 



Ai Iok 2 a iu)t 



2 £$»(» + D 

j n (ka) \h^_ x {kr)-^- r h^{kr) 

ht\\ka) j n -i(kr) - ^- j n (kr) 
I kr 



(50) 



and B v = 0, where the upper line holds for r > a and 
the lower line for r < a. As noted above, the coefficient 
a a in Eq. I|44|l indeed vanishes. Also, the discontinuity 
of the nth component of Bg in Eq. (|50(l at r — a is con- 
nected, according to Eq. with the nth component of 
the surface current density Js v obtained from Eqs. H48|) 
and l|33l) . A characteristic difference between this time- 
varying problem and the corresponding static case is the 
appearance of the spherical Bessel functions, which are 
solutions of the radial part of the Helmholtz equation in 
spherical coordinates. Using their limiting values [Tol ]. it 
can be seen that for k — > the static results in Eqs. i|35|) 
and (|36ll are obtained, as mathematically and physically 
expected. On the other hand, the radiative part of the 
external magnetic induction field, which decreases as 1/r, 
is given by 



(4n + 3)(2ra- 1)! 



B(Y 9 t) = 9 f J -° Ioka e i(kr-ua) 
l ' ' ' 4r 22«n!(n + l)! 

n=0 v ; 

xj 2n+1 (ka) Pl l+1 {cos 9). (51) 

In the dipole approximation, to < 1, this becomes the 
radiative magnetic induction field from an oscillating 
magnetic dipole of magnetic moment m = i\a 2 l i: 



B(r,t) 



Hok 2 

47T 



(fxm)xf 



i(kr— cjt) 



(52) 



The magnetic induction field in Eqs. (|49|l and (|50|l can be 

seen to be just that which is obtained with the more ar- 
duous technique of a dyadic Green's function expanded in 
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vector spherical harmonics and applied to the vector po- 
tential, which, by symmetry, only has the (^-component 
different from zero Q. As we have shown, a direct 
calculation of the electromagnetic field with r- and 9- 
components is much simplified if separation of variables 
is used. 



4. CONCLUSION 

For spherical coordinate systems, the Maxwell equa- 
tions outside sources lead to coupled equations involving 
all three components of the electromagnetic fields. In 
general, the statement is that one cannot separate spher- 
ical components of the Maxwell equations, and extensive 
techniques for solving the vector equations have been de- 
veloped which introduce vector spherical harmonics and 
use dyadic methods. We have shown, however, that sepa- 
ration of variables is still possible in the case of azimuthal 
symmetry, and so general solutions for each component of 
the electromagnetic vector fields were obtained. We have 



illustrated the use of these formulas with direct calcula- 
tions of electric and magnetic induction fields from local- 
ized charge and current distributions, without involving 
the electromagnetic potentials. Boundary conditions are 
easier to apply to these solutions, and their forms high- 
light the similarities and differences between the electric 
and magnetic cases in both time-independent and time- 
dependent situations. Finally, we remark that in cylin- 
drical coordinates, the other commonly used curvilinear 
coordinate system, the Maxwell equations do separate 
into equations for each vector component alone if there 
is cylindrical symmetry, so that the method of separation 
of variables can be used directly. 
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